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I.  RESULTS  FROM  PRIOR  SUPPORT 


1.  Simulation  Study  of  a  New  Mechanism  for  Excitation  of  Kinetic  Waves  in  a 
Magnetoplasma 


Nishikawa  et  al.  [1988]  have  investigated  the  new  ion-cyclotron-like  waves  bv  a  localized 
transverse  electric  field  by  means  of  simulation  with  the  assistance  of  the  nonlocal  kinetic 
theory.  The  linear  theory  shows  that  the  growth  rates  of  the  kinetic  Kelvin-Helmholtz  (K-H) 
modes  are  strongly  reduced  with  increasing  u  =  k^/ky,  and  they  become  unstable  only  when 
b  =  k2p2  <  1  and  kyL  s;  1,  where  L  is  the  scale  length  assocb' ->d  with  the  transverse  electric 
field.  On  the  other  hand,  the  new  modes  have  larger  frequencies  and  become  unstable  at  larger 
b  >  1  and  kyL  >  1  [Ganguli  et  al.,  1988]. 

Simulation  results  show  that  ion-cyclotron-like  waves  are  excited  in  regions  where  E  x  B 
drift  is  localized.  Linear  growth  rates  of  several  modes  are  obtained  from  the  wave  analysis 
of  the  simulation.  This  linear  analysis  shows  that  the  (0,  4)  mode  corresponds  to  large  6,  and 
large  kyL  has  the  maximum  growth  rate.  Clearly,  these  are  not  K-H  modes.  Further,  the 
simulation  results  show  that  density  gradients  help  to  enhance  the  growth  rates.  However,  like 
the  K-H  mode,  the  real  frequencies  of  this  instability  are  approximately  proportional  to  kyt  it- 
where  Vg  is  the  peak  value  of  the  E  X  B  drift. 

Nonlinear  phenomena  such  as  diffusion  and  coalescence  of  vortices  are  investigated.  In  the 
linear  stage  smaller  vortices  are  generated  and  larger  vortices  with  the  lower  frequencies  are 
dominant  in  the  nonlinear  stage.  In  the  nonlinear  stage  ions  diffuse  strongly  due  to  large-scale 
vortices. 

Recently,  W'e  have  investigated  the  electrostatic  waves  driven  by  the  combined  effects  of  a 
localized  transverse  electric  field  and  parallel  electron  drifts  by  means  of  simulation  with  the 
assistance  of  the  nonlocal  kinetic  theory  [Nishikawa  et  al.,  1989b, c,d]. 

We  have  performed  a  number  of  simulations  for  this  instability.  Simulation  results  show 
that  electrostatic  waves  are  excited  in  the  regions  where  the  E  X  B  drift  is  localized  in  the 
simultaneous  presence  of  parallel  electron  drifts  and  transverse  electric  fields.  Simulations  with 
only  parallel  electron  drifts  or  transverse  electric  fields  show  that  no  instability  grows  out  of 
the  thermal  noise.  The  simulation  with  both  the  parallel  electron  drift  and  the  transverse 
electric  field  shows  the  growing  waves  out  of  the  thermal  noise.  The  Doppler  shift  due  to  the 
E  x  B  drift  can  lower  the  phase  velocity  of  waves  along  the  magnetic  field.  Then  this  makes 
wave-particle  resonance  possible  for  smaller  V£,  which  leads  to  this  instability. 

The  nonlinear  phenomena  such  as  diffusion  and  coalescence  of  vortices  are  investigated 
[Nishikawa  et  al.,  1990a].  In  the  linear  stage,  smaller  vortices  are  generated  and  larger  vortices 
with  the  lower  frequencies  are  dominant  in  the  nonlinear  stage.  In  the  nonlinear  stage  the  ions 
diffuse  strongly  due  to  the  large  scale  vortices. 

The  spectrum  for  nonlinear  stage  was  analyzed  and  the  result  is  in  good  agreement  with 
the  satellite  data.  The  details  of  this  study  will  be  reported  in  J.  Geophys.  Res.  in  near  future. 
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2.  Beam  Instability  in  the  Foreshock 

As  an  application  of  the  simulation  method  used  in  the  proposed  research  (Broadband 
electrostatic  noise),  the  beam  instability  in  the  foreshock  has  been  investigated.  Electrons 
backstreaming  into  the  Earth’s  foreshock  generate  waves  near  the  plasma  frequency  by  the 
beam  instability.  Two  versions  of  the  beam  instability  exist:  the  ‘reactive’  version  in  which 
narrowband  waves  grow  by  bunching  the  electrons  in  space,  and  the  ‘kinetic’  version  in  which 
broadband  growth  occurs  by  a  maser  mechanism  [Cairns,  1887a,  b,  and  references  therein].  Re¬ 
cently,  Cairns  [1987b]  has  suggested  that  (1)  the  backstreaming  electrons  have  steep-sided  ‘cut¬ 
off’  distributions  which  arc  initially  unstable  to  the  reactive  instability,  (2)  the  back-reaction 
to  the  wave  growth  causes  the  instability  to  pass  into  its  kinetic  phase,  and  (3)  the  kinetic 
instability  saturates  by  quasilinear  relaxation. 

Cairns  and  Nishikawa  [1989]  have  performed  two-dimensional  simulations  of  the  reactive 
instability  for  Maxwellian  beams  and  cutoff  distributions.  The  results  of  the  simulations  are 
consistent  with  suggestions  (1)  and  (2)  above.  In  addition,  we  have  demonstrated  that  the 
reactive  instability  is  a  bunching  instability,  and  the  reactive  instability  saturates  and  passes 
over  into  the  kinetic  phase  bv  particle  trapping.  We  found  that  the  kinetic  growth  occurring 
after  saturation  of  the  reactive  instability  is  presumably  due  to  the  spatially  localized  gradients 
in  y  -  i’i;  base  space.  Both  simulation  results  and  numerical  solutions  of  the  dispersion 
equation  indicate  that  the  center  frequency  of  the  intense  narrowband  waves  near  the  foreshock 
boundary  may  be  between  0.9u;pe  and  0.98wpe,  rather  than  being  above  u  pe  as  previously 
believed.  Nonlinear  evolution  of  the  electron  plasma  waves  are  investigated  [Nishikawa  et  al. . 
1990c].  They  found  the  decay  process  in  which  a  large-amplitude  electron  plasma  wave  driven 
by  a  electron  beam  decays  into  a  backward  propagating  daughter  wave  and  an  ion  acoustic 
wave. 

Recent  one-dimensional  simulations  of  the  electron  beam-plasma  instability  [Dum,  1990] 
with  parameters  appropriate  for  the  electron  foreshock  quantitatively  tested  the  applicability 
of  quasi-linear  theory  and  other  turbulence  theories.  A  key  requirement  was  that,  in  contrast 
to  earlier  simulations,  the  model  at  all  times  allowed  the  excitation  of  a  large  number  of  wave 
modes.  Even  in  the  case  of  an  initially  cold  beam,  quasi-linear  theory  is  then  applicable, 
except  during  a  short  period  of  strong  beam  modulation  at  the  end  of  the  reactive  phase. 
Mode  coupling  is  seen  late  in  the  kinetic  phase  but  cannot  prevent  plateau  formation. 

It  is  important  to  investigate  possible  modifications  of  this  evolution  in  a  multidimensional 
system,  especially  for  a  low  temperature  beam,  which  initially  allows  the  excitation  of  a  large 
range  of  perpendicular  wave  numbers.  We  thus  compare  one-and  two-dimensional  simulations 
with  large  system  size.  The  effects  of  increased  magnetic  field  strength  (e.g.  along  auroral  field 
lines)  is  also  investigated.  This  study  will  be  reported  in  J.  Geophys.  Res. 
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3.  Whistler  Mode  Driven  by  the  Spacelab  2  Electron  Beam 

During  the  Spacelab  2  mission  while  an  electron  beam  was  being  ejected  from  the  shuttle, 
the  Plasma  Diagnostics  Package  (PDP)  detected  a  clear  funnel-shaped  emission  that  is  believed 
to  be  caused  by  whistler-mode  emission  from  the  electron  beam  [Gurnett  et  al. ,  1986].  In  order 
to  understand  the  mechanism  of  this  emission,  the  simulations  have  been  performed  using 
a  three-dimensional  magnetostatic  code  for  low-/?  plasmas  in  which  the  beam  electrons  arc 
initially  located  in  the  column  [Nishikawa  et  al.,  1989a],  In  order  to  simulate  the  continuous 
electron  ejection  from  the  shuttle,  the  simulations  were  also  performed  with  the  recycling  of 
the  beam  electrons.  The  beam  electrons  excite  whistler  waves  and  lower  hybrid  waves.  The 
brief  fluid  theory  based  on  the  magnetostatic  code  was  checked  with  the  simulation  results. 
The  propagating  whistler  mode  was  identified  with  the  theory.  The  simulation  results  show 
that  the  quasi  perpendicularly  (the  angle  between  the  magnetic  field  and  the  wavenumber  is 
larger  Ilian  50°)  propagating  whistler  waves  have  larger  amplitude  whose  real  frequencies  are 
smaller  than  the  local  electron  cyclotron  frequency.  This  fact  is  consistent  with  the  fact  that 
the  funnel-shaped  emission  is  observed  below  the  electron  cyclotron  frequency  away  from  the 
beam  electron.  The  beam  electrons  initially  in  the  column  diffuse  radially  as  well  as  slow  down 
due  to  the  E  X  B  caused  by  the  excited  waves.  The  acceleration  of  the  beam  electrons  also 
takes  place  due  to  the  excited  whistler  waves. 

In  order  to  compare  with  the  PDP  data,  the  local  magnetic  fields  Bx>y  and  the  perturbed 
electric  fields  EVtZ  are  diagnosed  at  the  several  points  in  the  simulation  system.  The  results 
show  that  the  waves  are  radially  excited  by  the  beam  electrons  localized  in  the  center  of  the 
system.  The  wave  spectra  of  the  electric  fields  Ey>z  diagnosed  at  x  =  31  Xe,y  ~  1 6 A e ,  and 
5  =  17Ae  show  the  several  kinds  of  waves  generated  by  the  beam  electrons.  The  analysis 
shows  that  the  lower  hybrid  waves  and  whistler  waves  are  excited.  The  frequency  range  of 
these  spectra  extends  from  uv,  to  beyond  u>pe  which  is  in  qualitative  agreement  with  the 
PDP  data.  As  the  PDP  data  show,  the  intense  electrostatic  narrowband  emission  around  the 
electron  plasma  frequency  has  been  observed  dominantly  in  the  spectra  of  the  electric  field  (  Ez ) 
with  and  without  the  recycling  of  the  beam  electrons.  This  wave  is  identified  as  the  parallel 
(quasi  parallel)  whistler  wave.  This  simulation  result  is  in  good  agreement  with  the  PDP  data 
[Gurnett  et  al.,  1986].  The  spectra  of  the  magnetic  fields  diagnosed  at  the  same  position  show 
that  whistler  modes  are  excited.  The  cB/E  was  calculated  for  the  whistler  mode.  In  the  case 
without  the  recycling  of  the  beam  electrons,  the  cBjE  is  approximately  0.35  for  the  whistler 
mode  (around  0.368wpe),  which  is  smaller  than  that  calculated  from  the  PDP  data  [Gurnett 
et  al.,  1986]. 

Recently  we  have  developed  a  3-D  electromagnetic  and  relativistic  particle  code  without 
solving  Poisson’s  equation.  “Radiating”  boundaries  are  introduced  besides  periodic  boundary 
conditions.  These  boundary  conditions  allow  radiation  to  escape  freely  out  of  the  simula¬ 
tion  domain  without  reflection.  Periodic  boundaries  are  applied  along  z  direction.  Radiating 
boundaries  are  used  for  the  boundaries  perpendicular  to  the  magnetic  field.  The  simulation 
domain  is  a  mesh  of  dimension  100  x  100  X  200  with  a  few  million  particles.  Whistler  emission 
as  shown  by  wave  data  obtained  from  the  Spacelab  2  mission  will  be  the  main  topics  of  this 
simulation  study. 
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4.  Simulation  of  Electron  Cyclotron  Harmonic  Waves  by  AMPTE  Lithium  Release 

We  have  studied  the  generating  electrostatic  turbulence  in  the  lithium  cloud  due  to  solar 
wind  particles  which  are  flowing  through  it  [Nishikawa  et  al.,  1990b].  Several  simulations  have 
been  performed  in  order  to  understand  the  processes  generating  electron  cyclotron  harmonic 
waves. 

The  simulations  show  that  electron  cyclotron  harmonic  waves  are  excited  by  mainly  the 
solar  wind  electrons.  The  spectra  of  perturbations  of  electrostatic  potential  exhibit  electron 
cyclotron  harmonic  waves. 

The  overall  frequency  range  is  obtained  by  observing  the  electric  fields  at  some  points  in 
the  simulation  region.  Power  spectra  of  the  electric  fields  exhibit  electron  cyclotron  frequency 
and  higher  harmonics  which  is  simular  to  the  spectrum  obtained  by  the  satellite. 

Possible  mechanisms  of  generating  electron  cyclotron  harmonic  waves  have  been  discussed 
in  the  previous  work  [Roeder  et  al.,  1987].  In  this  simulation  study,  the  possible  energy 
sources  of  the  instability  are  included  automatically.  The  relative  velocities  among  different 
species  become  the  energy  source  of  instabilities.  The  relative  velocity  between  solar  wind 
electrons  and  lithium  ions  would  be  one  possible  energy  source  for  the  excitation  of  ECU 
waves.  The  fact  that  the  instability  can  be  excited  without  new-born  cold  electrons  support 
this  mechanism.  Theoretically,  the  mode  couplings  among  harmonics  occur  »  ainly  due  to  the 
electron  perpendicular  drift,  which  lead  to  the  instability  [Tataronis  and  Crawford,  1970]. 

According  to  Figure  1  by  Gurnett  et  al.  [1986b],  the  higher  frequency  ranges  of  waves  look 
very  simular  at  both  cases  with  as  npr  and  %  <  npr.  In  this  study,  we  found  that  ECU 
waves  are  excited  also  in  the  case  of  nu  =  npr.  The  satellite  data  show  that  ECU  waves  can  be 
observed  only  in  some  range  of  lithium  ion  density  (0.03  -  0.18  cm-3)  associated  with  lithium 
drift  energy  [Roeder  et  al.,  1987].  The  observations  show  that  in  the  case  of  %  ss  npr  ECU 
waves  are  not  found.  The  simulation  results  show  that  in  the  same  case  of  riu  =  npr  ECU  waves 
are  excited.  This  discrepancy  may  be  explained  by  the  setting  of  simulation  geometry  due  to 
the  two-dimensional  code.  If  the  simulation  plane  is  set  nearly  parallel  to  the  magnetic  fields 
(Boy/ B0  <  1),  the  angle  between  the  magnetic  fields  and  the  drift  velocity  of  hot  electrons  and 
protons  becomes  very  small.  Then,  in  the  case  of  n;,  =  npr,  the  parallel  drift  of  solar  wind 
protons  and  electrons  may  excite  ion  acoustic  waves  more  dominantly  than  ECU  waves,  which 
is  not  included  in  this  study  [Gurnett  et  al.,  1986b;  Ma  et  al.,  1987]. 

Further  study  is  required  on  the  observed  plasma  velocity  distributions  and  comparisons 
of  simulation  wave  results  with  the  wave  power  spectra.  More  simulations  are  required  for 
further  progress  in  the  interpretation  of  electrostatic  noise. 
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